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Emergence  and  Frustration  of 
Magnetism  with  Variable-Range 
Interactions  in  a  Quantum  Simulator 

R.  Islam,1*  C.  Senko,1  W.  C.  Campbell,1!  S.  Korenblit,1  ].  Smith,1  A.  Lee,1  E.  E.  Edwards,1 
C.-C.  ].  Wang,2!  J-  K.  Freericks,2  C.  Monroe1! 

Frustration,  or  the  competition  between  interacting  components  of  a  network,  is  often  responsible 
for  the  emergent  complexity  of  many-body  systems.  For  instance,  frustrated  magnetism  is  a 
hallmark  of  poorly  understood  systems  such  as  quantum  spin  liquids,  spin  glasses,  and  spin 
ices,  whose  ground  states  can  be  massively  degenerate  and  carry  high  degrees  of  quantum 
entanglement.  Here,  we  engineer  frustrated  antiferromagnetic  interactions  between  spins  stored 
in  a  crystal  of  up  to  16  trapped  171Yb+  atoms.  We  control  the  amount  of  frustration  by 
continuously  tuning  the  range  of  interaction  and  directly  measure  spin  correlation  functions 
and  their  coherent  dynamics.  This  prototypical  quantum  simulation  points  the  way  toward  a  new 
probe  of  frustrated  quantum  magnetism  and  perhaps  the  design  of  new  quantum  materials. 


Predicting  the  behavior  of  many-body  quan¬ 
tum  materials  such  as  frustrated  magnets 
is  difficult  because  the  number  of  relevant 
configurations  scales  exponentially  with  the  num¬ 
ber  of  particles  (1-3).  Feynman  proposed  the  use 
of  a  quantum  simulator  for  the  task.  Here,  in¬ 
teractions  are  engineered  in  a  “standard”  quan¬ 
tum  system  to  illuminate  the  physics  behind  the 
real  material  (4).  Cold  atoms  are  excellent  stan¬ 
dards  for  quantum  simulations  of  magnetism, 
with  the  ability  to  tailor  frustrated  magnetic  in¬ 
teractions  and  measure  the  individual  atomic 
spins  (5,  6).  Neutral  atomic  systems  are  typi¬ 
cally  limited  to  nearest-neighbor  interactions  (7), 
although  geometrically  frustrated  interactions  can 
be  realized  in  certain  optical  lattice  geometries 
(8).  The  strong  Coulomb  interaction  between  cold 
atomic  ions  (9)  has  led  to  the  realization  of  long- 
range  Ising  couplings  between  individual  trapped 
ion  spins  (10-14)  and  the  observation  of  spin  frus¬ 
tration  and  quantum  entanglement  in  the  smallest 
system  of  three  spins  (75). 

Here,  we  report  the  implementation  of  variable- 
range  antiferromagnetic  (AFM)  Ising  interactions 
with  up  to  16  atomic  ion  spins,  using  optical  di¬ 
pole  forces.  We  directly  measure  the  emergence 
and  frustration  of  magnetic  ordering  through  spa¬ 
tially  resolved  imaging  of  the  ions.  The  spins  are 
initially  polarized  along  a  strong  effective  mag¬ 
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netic  field  transverse  to  the  Ising  couplings,  and 
when  the  field  is  lowered,  the  spins  order  them¬ 
selves  according  to  the  characterstics  of  the  Ising 
interactions.  We  can  increase  the  level  of  frustra¬ 
tion  by  increasing  the  range  of  the  interactions, 
which  results  in  a  more  equitable  balance  of  com¬ 
peting  interactions  and  a  suppression  of  magnetic 
order.  The  quantum  coherence  in  the  system  is 
characterized  by  reversing  the  transverse  field 
back  to  its  initial  value  and  comparing  the  result¬ 
ing  state  with  the  initial  state.  These  experiments 
present  simulations  in  a  nontrivial  quantum  sys¬ 
tem  that  approaches  a  complexity  level  at  which 
it  becomes  difficult  or  impossible  to  calculate  the 
spin  dynamics. 


Model 

We  simulate  the  quantum  transverse  Ising  mod¬ 
el  with  long-range  AFM  interactions  in  a  one¬ 
dimensional  spin  chain,  given  by  the  Hamiltonian 


where  the  Planck  constant  h  is  set  to  1 ,  (y  = 

x,  y,  z)  are  the  spin- 1/2  Pauli  operators  for  the 
z'th  spin  (i  =  1,2,... A),  B  is  the  effective  trans¬ 
verse  magnetic  field,  and  Jy  >  0  is  the  Ising  cou¬ 
pling  between  spins  i  and  y,  falling  off  with  the 
lattice  spacing  |  i  -j\  approximately  as 


where  0  <  a  <  3  (5). 

For  B  »  Jy  on  all  pairs,  the  spins  are  po¬ 
larized  along  the  effective  transverse  magnetic 
field  in  the  ground  state  \^y  ]y  ]y  ...)  of  the 
Hamiltonian  in  Eq.  1 ,  where  |  ]y)  denotes  a  spin 
along  the  +y  direction  of  the  Bloch  sphere.  As 
the  ratio  of  B  to  the  Ising  couplings  is  reduced, 
the  system  crosses  over  to  an  ordered  state  dic¬ 
tated  by  the  form  of  the  Ising  couplings,  and  the 
spectrum  of  energy  levels  depends  on  the  range 
of  the  interactions.  For  any  finite  a  >  0,  we  find 
by  direct  calculation  that  the  staggered  AFM  states 
\UU  ...)  and  lltlt  •••)  constitute  the  doubly 
degenerate  ground  state  manifold  at  B  =  0,  with 
the  degeneracy  arising  from  the  time  reversal  or 
the  global  spin  flip  symmetry  of  the  Hamiltonian. 
Here  ||)  and  ||)  are  spin  states  oriented  along  the 
Ising  or  x  direction  of  the  Bloch  sphere.  Thus,  the 
system  exhibits  nearest-neighbor  AFM  or  Neel 
ordering  at  sufficiently  low  temperatures.  When 
the  interactions  are  uniform  over  all  pairs  of  spins 
(a  — >  0),  the  system  becomes  maximally  frustrated 
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Fig.  1.  Theoretical  energy  spectrum  and  critical  gap  in  the  long-range  antiferromagnetic  Ising 
model  (Eq.  1)  for  N  =  10  spins.  (A)  Low-lying  energy  states  for  a  =  1  (characteristic  range  of  £  =  5  sites) 
as  a  function  of  the  dimensionless  parameter  Bl]0.  The  spacing  between  the  ground  state  at  E  =  E0  and  the 
first  coupled  excited  state  (black  lines)  reaches  a  bottleneck  at  a  critical  value  Bcl]0  with  critical  gap  Ac.  (B) 
Theoretical  dependence  of  Bcl]0  (dotted  line)  and  Acl]0  (solid  line)  on  the  range  of  the  interaction.  As  the 
interaction  range  increases,  the  competing  long-range  couplings  make  it  easier  to  create  excitations  and 
the  critical  gap  is  reduced,  so  a  relatively  small  effective  transverse  field  can  break  the  spin  ordering.  Both 
parameters  approach  zero  as  a  -►  0  or  £  -►  oo.  Present  experiments  are  performed  with  parameters  in 
the  shaded  region. 
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and  the  excitation  gap  (Fig.  1  A)  closes,  leading 
to  a  finite  entropy  density  in  the  ground  state. 
In  this  case,  any  spin  configuration  with  a  net 
magnetization  of  zero  (1/2)  for  even  (odd)  num¬ 
bers  of  spins  belongs  to  the  ground  state. 

Between  the  limits  B  =  0  and  B  »  Jy,  the 
energy  spectrum  features  a  minimum  gap,  whose 
position  and  size  depends  on  the  degree  of  frus¬ 
tration  in  the  system,  or  the  interaction  range. 
(The  interaction  range  is  defined  as  the  number  of 
lattice  spacings  \  where  the  interaction  falls  off  to 
20%  of  the  nearest-neighbor  Ising  coupling:  \  = 
51/ct.)  Figure  1A  shows  a  few  low-lying  energy 
states  of  the  Hamiltonian  in  Eq.  1  for  an  in¬ 


teraction  range  of  \  =  5  (corresponding  to  a  =  1). 
The  first  excited  eigenstate  merges  with  the  ground 
state  for  small  B/J0  and  has  the  same  spin  order  as 
the  ground  state  near  B/J0  =  0.  The  critical  gap  Ac 
between  the  ground  and  the  first  coupled  excited 
state  determines  the  adiabaticity  criterion  (16). 
Figure  IB  compares  the  position  (dotted  line) 
and  size  (solid  line)  of  the  critical  gap  of  the 
Hamiltonian  for  various  ranges.  As  the  range  and 
hence  the  amount  of  frustration  increases,  the  crit¬ 
ical  field  is  pushed  toward  zero,  and  the  gap 
closes.  To  observe  the  effects  of  frustration,  re¬ 
flected  in  the  density  of  states  near  the  ground 
state,  we  quench  the  system  by  ramping  the  ef¬ 


fective  transverse  magnetic  field  faster  than  the 
critical  gap  (\B/B\  >  Ac)  to  populate  the  lowest 
coupled  excited  states.  The  observed  spin  order 
depends  on  the  resulting  degree  of  excitation  and 
hence  on  the  level  of  frustration. 

Experiment 

The  spins  are  realized  in  a  collection  of  171Yb+ 
ions  confined  in  a  linear  radiofrequency  (Paul) 
trap,  with  the  effective  spin- 1/2  system  repre¬ 
sented  by  two  hyperfme  “clock”  states  within 
each  ion  |fz)  and  |  jz),  separated  by  the  hyper¬ 
fme  frequency  vHf  =  12.642819  GHz  (1 7).  The 
variable-range  AFM  Ising  interactions  are  gen¬ 
erated  by  applying  off-resonant  spin-dependent 
optical  dipole  forces  (11)  that  drive  stimulated 
Raman  transitions  between  the  spin  states  while 
modulating  the  Coulomb  interaction  between  the 
ions  in  a  controlled  way  (18).  The  effective  mag¬ 
netic  field  is  generated  by  simultaneously  driving 
coherent  transitions  between  the  spin  states  with 
a  7t/2-phase  shift  with  respect  to  the  dipole  force 
beams.  At  any  time,  we  measure  the  state  of  the 
spins  by  illuminating  the  ions  with  resonant  ra¬ 
diation  and  collecting  state-dependent  fluorescence 
on  an  imager  with  site-resolving  optics  (1 7).  From 
this  information,  we  can  extract  all  spin  correla¬ 
tion  functions  (18). 

The  quantum  simulation  begins  by  optically 
pumping  all  spins  to  the  |  jz)  state  and  then  co¬ 
herently  rotating  them  all  about  the  x  axis  of  the 
Bloch  sphere  to  initialize  each  spin  in  state  \\y) 
along  the  effective  transverse  magnetic  field. 


Fig.  2.  Spin  order  versus 
speed  of  ramp,  for  l\l=  10 
spins.  The  spins  are  initial¬ 
ized  with  B/J o  =  5  and  the 
transverse  field  is  ramped 
exponentially  down  for  six 
time  constants,  and  the  ex¬ 
periment  is  repeated  for  var¬ 
ious  values  of  time  constant  x. 

Symbols  (solid  line)  indicate 
the  scaled  staggered  Binder 
cumulant  gSf  versus  the  to¬ 
tal  duration  6x  of  the  ramp  _  .  _  A  _ 

measured  (expected  from  Duration  of  quantum  simulation  (=6x,  ms) 

theory).  As  the  Hamiltonian  evolves  more  slowly,  the  observed  spin  order  shows  more  ground  state 
order,  and  less  excitation  for  ramp  times  under  -  2.5  ms.  For  longer  times,  the  spins  become  disordered, 
implying  external  decoherence  in  the  system.  Here  (and  in  the  following  figures)  the  error  bars  include 
statistical  fluctuations  and  estimated  detection  uncertainties. 


Fig.  3.  Quantum  phase  transition  from  a  paramagnet  to  an  antifer- 
romagnet  (Eq.  1),  with  ]y  «  Ir  -  j I-1*05  in  a  system  of  10  spins.  The 

exponent  (a  =  1.05)  is  estimated  from  the  average  couplings  between  spins  in 
this  inhomogeneous  system.  (A)  Image  of  10  trapped  ions,  with  a  distance  of 
22  pm  between  the  first  and  last  ion.  (B)  Measured  two-point  correlation 
function  between  a  chosen  spin  (on  the  edge)  and  the  other  spins  separated 
by  r  lattice  sites,  C1/1+r  =  (a^o£1+r))  -  (cri1^1^),  averaged  over  4000  ex¬ 
periments  for  each  value  of  the  parameter  B/J0.  For  B/J0  =  5,  the  spins  are 
initially  polarized  along  the  transverse/ field  with  little  correlation  along 
the  Ising  x  direction.  As  the  field  is  reduced,  the  spins  cross  over  to  pre¬ 


dominantly  AFM  states  I T  i  t  i  . . .)  and  U  t  i  t  . . .),  resulting  in  alternating 
signs  in  the  two-point  correlations  with  separation.  The  solid  lines  are  shown  to 
guide  the  eye.  (C)  Measured  occurrence  probability  of  all  210  =  1024  states  at 
B/J0  «  5  (paramagnetic  state,  red  trace)  and  B/J0  «  0.01  (AFM  phase,  black 
trace).  The  states  are  listed  in  binary  order,  with  spin  U)  =  0  and  1 t )  =  1.  The 
residual  peaks  in  the  red  trace  are  consistent  with  detection  errors  biased 
toward  states  with  many  It)  spins  such  as  127,  511,  and  1023.  The  two  tall 
peaks  in  the  black  trace  at  341  (0101010101)  and  682  (1010101010) 
correspond  to  two  Neel-ordered  staggered  AFM  states,  shown  with  camera 
images  of  these  cases  and  contributing  -17%  to  the  population. 
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The  Hamiltonian  (Eq.  1)  is  then  switched  on  with 
an  initial  field  B0  ~  5  J0 ,  where  J0  (~1  kHz)  is 
the  average  nearest-neighbor  Ising  coupling,  thus 
preparing  the  spins  in  the  ground  state  of  the 
initial  Hamiltonian  with  a  fidelity  better  than 
97%.  The  effective  magnetic  field  is  ramped 
down  exponentially  in  time  with  a  time  constant 
of 400  ps  to  a  final  value  B  of  the  transverse  field, 
but  no  longer  than  2.4  ms  overall,  to  avoid  de¬ 
coherence  effects  present  in  the  system  (see 
Fig.  2).  At  this  point,  the  Hamiltonian  is  switched 
off,  freezing  the  spins  for  measurement.  We  then 
measure  the  x  or  y  component  of  each  spin  o® 
or  oW  by  first  rotating  our  measurement  axes 
with  an  appropriate  global  te/2  pulse  similar  to 
the  initialization  procedure,  before  capturing 
the  spin-dependent  fluorescence  on  the  im¬ 
ager.  The  experiments  are  repeated  2000  to  4000 
times  to  measure  expectation  values  of  certain 
spin  operators  and  correlation  functions  (denoted 
by  <... )). 

Order  Parameters  and  Correlation  Functions 

From  these  measurements,  we  can  construct 
order  parameters  appropriate  for  observing  low- 
energy  AFM  states.  Various  moments  of  the 
staggered  magnetization  operator 

N 

ms  =  2  (-l)zo^  distinguish  between  para¬ 

magnetic  and  AFM  order,  and  also  quantify  spin 
flip  excitations.  In  particular,  the  normalized  fourth 
moment  of  this  magnetization,  known  as  the 
Binder  cumulant  (19)  g,  =  ((mg  -  ( m,))4)/((m v  - 
<ms»2)2,  scaled  togs  =  (g°s  -  gs)/(g°  ~  1)  to 
remove  finite  size  effects,  varies  from  gs  =  0 
to  gs  =  1  as  the  paramagnetic  state  gives  way 
to  AFM  order  (see  Fig.  2).  Here  g®  =  3  -  2/ N 
is  the  Binder  cumulant  in  a  perfect  paramag¬ 
netic  state  of  the  N  spins.  We  can  also  form  any 
correlation  function  of  the  spins  such  as  the  two- 
point  correlation  Qj  =  (a$  a® )  -(a^)(a^),  al¬ 
lowing  a  direct  probe  of  spin  order  for  each 
experimental  realization.  The  Fourier  transform 


of  this  correlation  function  is  the  structure  func 
tion  S(k)  =  |  X  C(r)elkr\  ,  where  C(r)= 

N-r  r=l 

2  Cm>w+r  is  the  average  correlation  over  r 

sites  in  the  chain.  The  structure  function  shows 
spin  order  versus  wave  number  k,  with  S(k  =  n) 
singling  out  the  presence  of  the  nearest-neighbor 
Neel  AFM  order. 

Figure  3  shows  the  onset  of  AFM  ordering  in 
the  quantum  simulation  of  the  frustrated  trans¬ 
verse  field  Ising  model  in  a  system  of  10  spins. 
Two-point  spin  correlations  C^r  between  a  cho¬ 
sen  edge  spin  and  the  other  spins  are  presented  in 
Fig.  3B  at  various  stages  in  the  ramp  B/J0  =  5, 
0.25,  and  0.01.  For  larger  transverse  magnetic 
fields,  there  are  no  appreciable  correlations  be¬ 
tween  the  spin  components  along  the  Ising  di¬ 
rection.  As  the  ratio  of  B/J0  is  lowered,  however, 
a  zig-zag  pattern  emerges,  with  negative  (posi¬ 
tive)  correlations  between  spins  separated  by  odd 
(even)  lattice  spacings.  For  B/J0  ~  0.01,  the 
nearest-neighbor  spin  correlation  reaches  about 
60%,  and  the  correlation  length  (defined  to  be  the 
distance  at  which  the  absolute  correlation  drops 
to  He  of  the  nearest  neighbor  value)  reaches 
about  six  lattice  sites.  The  effective  field  was 
ramped  exponentially  down  from  B/J0  ~  5  with  a 
time  constant  of  400  ps  in  this  experiment.  This 
ramping  is  not  slow  enough  to  be  adiabatic,  and 
the  diabatic  effects  prevent  the  spin  ordering  from 
reaching  a  perfect  AFM  phase.  Figure  3C  shows 
the  measured  probabilities  of  all  2 10  =  1024  pos¬ 
sible  spin  states  measured  along  the  Ising  x  direc¬ 
tion,  sorted  in  binary  order  with  spin  |j)  =  0  and 
It)  =  1.  The  net  detection  fidelity  of  each  spin  is 
-97%,  after  postfiltering  the  measurements  based 
on  calibrating  the  known  detection  errors  for  each 
spin  (18, 20).  The  initial  paramagnetic  phase  shown 
in  red  (BIJq  ~  5)  exhibits  a  roughly  uniform 
probability  of  1/1024  ~  0.1%  for  each  state  (the 
residual  peaks  in  the  red  trace  are  consistent  with 
detection  errors).  The  spin-ordered  phase  shown 


in  black  (BIJq  =  0.01)  displays  the  emergence 
of  the  two  AFM  states,  each  with  an  occupation 
probability  of  about  8.5%.  Other  prominent  peaks 
correspond  to  single  spin-flips  and  other  low- 
lying  excitations  from  the  two  ground  states. 

In  Fig.  4,  we  probe  the  frustration  in  the  sys¬ 
tem  for  various  ranges  of  interactions.  Here,  we 
look  at  the  spin  order  achieved  in  the  quantum 
simulation  when  the  external  magnetic  field  is 
ramped  down  to  B/J0  ~  0.01  for  four  different 
ranges  of  interactions.  In  Fig.  4A,  we  show  the 
measured  structure  function  S(k)  at  wave  vector 
values  k  =  ,  fg , ....  tt  from  the  measured  two- 

point  correlation  functions.  To  directly  compare 
the  different  interaction  ranges,  we  choose  the 
same  external  magnetic  field  ramp  time  constant, 
t  =  0.4/Jq-  As  the  range  of  interaction  increases, 
the  ground  state  AFM  order  (given  by  the  struc¬ 
ture  function  peak  at  k  =  n)  disappears,  reflecting 
increased  occupation  of  the  excited  states  as  the 
frustration  grows.  Figure  4B  displays  the  ob¬ 
served  distribution  of  energy  P(Ei)  for  a  =  1.05 
(shorter  range)  and  a  =  0.76  (longer  range)  power- 
law  exponents,  along  with  the  cumulative  energy 
distribution  function.  The  eigenenergies  of  each 
configuration  are  calculated  using  Eq.  1  with 
B  =  0.  For  the  longer-range  interactions,  exci¬ 
tations  are  more  prevalent,  and  the  energy  gap 
between  the  ground  and  the  first  excited  state  is 
reduced,  both  of  which  are  signatures  of  increas¬ 
ing  frustration  in  the  system.  The  observed  final 
entropy  per  particle  S  =  ~^P(Ei)\ogP(Ei)  is 
seen  to  increase  from  0.832  to  0.903  as  the  in¬ 
teraction  range  grows  from  a  =  1 .05  (£,  =  4.6  sites) 
to  a  =  0.67  (£  =  11  sites),  which  is  also  a  signa¬ 
ture  of  the  increased  frustration  in  the  system.  As 
a  reference,  the  paramagnetic  state  distribution 
shows  an  entropy  per  particle  of  0.959,  which  is 
slightly  less  than  unity  because  of  detection  errors. 

Quantum  Coherence 

The  above  measurements  of  the  state  distribu¬ 
tion  concern  only  the  diagonal  components  or 
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Fig.  4.  (A)  Structure  function  S(k)  versus  wave  vector  k  for  various  ranges 
of  AFM  interactions,  for  B/]0  =  0.01  in  a  system  of  N  =  10  spins.  The 

increased  level  of  frustration  for  the  longer-range  interactions  reduces  the 
observed  antiferromagnetic  spin  order.  The  detection  errors  may  be  larger  than 
shown  here  for  the  longest  range  of  interactions,  owing  to  spatial  crosstalk 


from  their  closer  spacing.  (B)  Distribution  of  observed  states  in  the  spin  sys¬ 
tem,  sorted  according  to  their  energy  £,  that  was  previously  calculated  by 
diagonalizing  Eq.  1  with  B  =  0.  Data  are  presented  for  two  ranges  (red  for 
a  =  1.05  and  blue  for  a  =  0.76).  The  dashed  lines  indicate  the  cumulative 
energy  distribution  functions  for  these  two  ranges. 
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populations  of  the  density  matrix.  To  character¬ 
ize  the  quantum  coherence  in  the  simulation,  we 
retrace  the  external  magnetic  field  back  to  its  ini¬ 
tial  value  after  ramping  it  down  to  almost  zero 
and  measure  each  spin  along  the  transverse  (y) 
magnetic  field.  Figure  5  shows  the  distribution  of 
the  measured  value  my  of  the  total  transverse  spin 
operator^  =  X  oj^ at  three  different  times:  first 
at  the  beginning  of  the  simulation,  second  when 
the  transverse  field  has  been  ramped  down  to 
nearly  zero,  and  third  after  the  transverse  field 
returns  to  its  initial  value  (a  =  0.92  for  these 
data).  The  initial  state  is  ideally  a  delta  function 
at  (Sy)  =10,  but  finite  detection  efficiency  broad¬ 
ens  the  distribution.  At  the  lowest  value  of  the 
field  ( B/J0  ~  0.01),  the  transverse  magnetiza¬ 
tion  is  distributed  near  (Sy)  =  0,  as  the  spins  are 
presumably  ordered  along  the  Ising  x  direction. 
When  the  external  field  is  ramped  back  to  its  ini¬ 
tial  state,  the  distribution  of  the  total  spin  returns 
toward  the  initial  distribution,  with  an  average  mag¬ 
netization  that  is  approximately  (Sy)  =  68(4)%  of 


its  initial  value,  indicating  that  some  level  of 
quantum  coherence  is  maintained  throughout  the 
simulation.  This  number  is  in  agreement  with  the 
theoretically  estimated  average  magnetization  of 
(Sy)  ~  70%  of  its  initial  value,  obtained  by  nu¬ 
merically  integrating  the  Schrodinger’s  equation 
without  any  decoherence. 

To  probe  decoherence  in  the  simulation,  we 
repeat  the  experiment  with  various  ramping  speeds 
of  the  effective  magnetic  field.  In  Fig.  2,  we  plot 
the  AFM  order  parameter  gs  versus  the  total  du¬ 
ration  for  the  experiment  for  a  long-range  cou¬ 
pling  (a  =  1.05)  for  N=  10  spins.  Each  data  point 
represents  the  spin  order  achieved  after  ramping 
the  magnetic  field  B  down  exponentially  from 
5Jq  for  a  total  duration  of  six  time  constants.  The 
AFM  order  grows  with  slower  ramping,  as  ex¬ 
pected,  for  up  to  t  —  400  jlus  .  But  we  also  observe 
a  saturation  and  then  decay  in  the  spin  order, 
which  might  indicate  the  presence  of  decoher¬ 
ence  in  the  system  at  long  times.  During  the  sim¬ 
ulation,  spontaneous  Raman  scattering  from  the 


optical  beams  is  expected  to  occur  at  a  rate  of  less 
than  6  s~]  per  spin  (27),  which  is  consistent  with 
separate  measurements  of  the  spin  relaxation 
from  a  single  spin  and  is  therefore  not  expected 
to  contribute  to  decoherence  given  the  time  scales 
in  the  experiment.  The  phonon  population  is 
expected  to  be  well  under  10%  for  all  the  data 
presented  here  (22).  A  principal  source  of  deco¬ 
herence  appears  to  be  the  intensity  fluctuations 
in  the  Raman  beams,  from  beam  pointing  in¬ 
stabilities  and  fluctuations  in  the  optical  power. 

Maintaining  adiabaticity  while  ramping  the 
magnetic  field  down  is  more  difficult  when  the 
Ising  interactions  in  the  experiment  are  frustrated, 
because  the  relevant  energy  gaps  are  smaller.  To 
directly  compare  frustrated  versus  nonffustrated 
systems,  we  execute  quantum  simulations  of  both 
long-range  AFM  and  ferromagnetic  (FM)  Ising 
models  in  a  system  of  N=  16  spins  (Fig.  6).  For 
the  FM  experiment  (Fig.  6C),  we  initialize  the 
spins  in  the  highest-energy  state  with  respect  to 
the  transverse  field  . . .)  and  ramp  the  field 

down  as  before  (10,  15).  For  the  AFM  experi¬ 
ment  with  the  same  ramp  rate,  we  find  that  the  best 
AFM  nearest-neighbor  correlation  (Fig.  6 A)  is 
only  30(3)%,  corresponding  to  a  staggered  mag¬ 
netization  of  about  30(2)%,  whereas  the  simu¬ 
lation  of  the  FM  model  shows  a  clear  FM  spin 
order  across  the  chain  (Fig.  6B),  reaching  73(3)% 
magnetization.  We  have  also  observed  a  level  of 
70(10)%  FM  magnetization  emerging  in  N=  18 
spins  (18). 

The  interacting  spin  system  that  we  study  is 
approaching  a  complexity  level  at  which  it  be¬ 
comes  difficult  or  impossible  to  predict  its  be¬ 
havior.  Static  properties  such  as  the  ground  state 
order  or  the  excited  state  energies  can  be  cal¬ 
culated  for  hundreds  of  spins  using  Monte  Carlo 
methods  (23);  however,  the  calculation  of  dynam¬ 
ics  of  fully  connected  spin  models  is  currently 
limited  to  approximately  N  =  30  spins  (24,  25). 
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Fig.  5.  Quantum  coher-  1 . 0 

ence  in  the  simulation, 
in  a  system  of  N  =  10 
spins.  Probabilities  of 
different  values  of  the 
total  spin  component 
along  y  in  the  initial 
polarized  state  (red), 
when  the  transverse  field 
is  ramped  to  «0.01  ]0 
(black),  and  when  the 
transverse  field  is  reversed 
back  to  its  initial  value 
(green).  After  reversal,  the 
y  magnetization  returns 

to  ~68(4)%  of  its  initial  value,  indicating  quantum  coherence  in  the  evolution.  The  trajectory 
transverse  field  ( B ,  in  green)  and  all  the  Ising  couplings  (],  in  blue)  are  shown  in  the  inset,  a  =  0 
these  data. 
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Fig.  6.  Magnetic  or¬ 
dering  in  N=  16  spins. 

(A)  Image  of  16  trapped 
ions,  with  a  distance  of 
30  pm  between  the  first 
and  last  ion.  (B)  Pair  cor¬ 
relation  function  mea¬ 
sured  at  various  stages 
of  the  quantum  simula¬ 
tion,  for  B/J0  =  5  (red) 
and  B/Jo  =  0.01  (blue) 
in  an  AFM  coupling  fall¬ 
ing  off  with  distance  as 
Jij  ~  1/  -  jT1  among  N  = 

16  spins.  Small  amounts 
of  staggered  order  are 
seen,  tempered  by  the 
small  gaps  and  frustra¬ 
tion  in  the  low-energy 
states.  (C)  In  contrast, 

for  all  FM  couplings  (again  with  77y  ~  1/  - ; I-1),  the  gaps  are  large  and  clear  FM 
order  is  seen.  Here  the  measured  distribution  of  magnetization  is  plotted.  The 
paramagnetic  phase  of  16  spins  is  indicated  in  red,  and  after  the  field  is 
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ramped  to  nearly  zero,  the  distribution  clearly  bifurcates,  indicating  popula¬ 
tion  weighted  heavily  toward  the  FM  states  1 1  i  i . . .)  and  1 1 1 1 . . .),  indicated 
in  blue.  The  resulting  magnitude  magnetization  is  -73%. 
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Small  ion  trap  quantum  simulators  such  as  that 
reported  here  may  soon  reach  this  milestone  with 
technical  upgrades  in  the  hardware,  including 
lower  vacuum  chamber  pressures  to  prevent  col¬ 
lisions  with  the  background  gas,  better  stability  of 
the  optical  intensities,  and  higher  optical  power 
so  that  fluctuations  in  the  beam  inhomogeneities 
can  be  suppressed. 
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Kepler-62:  A  Five-Planet  System 
with  Planets  of  1.4  and  1.6  Earth 
Radii  in  the  Habitable  Zone 
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Thomas  N.  Gautier  III,22  John  C.  Geary,3  Alan  Gould,23  Michael  Haas,1  Christopher  Henze,1 
Andrew  W.  Howard,24  Steve  B.  Howell,1  Daniel  Huber,1  Jon  M.  Jenkins,5  Hans  Kjeldsen,18 
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We  present  the  detection  of  five  planets — Kepler-62b,  c,  d,  e,  and  f — of  size  1.31,  0.54,  1.95, 
1.61  and  1.41  Earth  radii  (/?0),  orbiting  a  K2V  star  at  periods  of  5.7,  12.4,  18.2,  122.4,  and 
267.3  days,  respectively.  The  outermost  planets,  Kepler-62e  and  -62f,  are  super-Earth-size 
(1.25  /?0  <  planet  radius  <  2.0  /?0)  planets  in  the  habitable  zone  of  their  host  star,  respectively 
receiving  1.2  ±  0.2  times  and  0.41  ±  0.05  times  the  solar  flux  at  Earth's  orbit.  Theoretical  models 
of  Kepler-62e  and  -62f  for  a  stellar  age  of  ~7  billion  years  suggest  that  both  planets  could 
be  solid,  either  with  a  rocky  composition  or  composed  of  mostly  solid  water  in  their  bulk. 


Kepler  is  a  NASA  Discovery-class  mission 
designed  to  determine  the  frequency  of 
Earth-radius  planets  in  and  near  the  hab¬ 
itable  zone  (HZ)  of  solar-like  stars  (1-6).  Planets 
are  detected  as  “transits”  that  cause  the  host  star  to 
appear  periodically  fainter  when  the  planets  pass 
in  front  of  it  along  the  observer’s  line  of  sight. 
Kepler-62  [Kepler  Input  Catalog  (KIC)  9002278, 
Kepler  Object  of  Interest  (KOI)  701]  is  one  of  about 


170,000  stars  observed  by  the  Kepler  spacecraft. 
On  the  basis  of  an  analysis  of  long-cadence  photo¬ 
metric  observations  from  Kepler  taken  in  quarters 
1  through  12  (13  May  2009  through  28  March 
2012),  we  report  the  detection  of  five  planets  or¬ 
biting  Kepler-62,  including  two  super-Earth-size 
planets  in  the  HZ  as  well  as  a  hot  Mars-size  planet 
(Fig.  1  and  Table  1).  Before  validation,  three  of 
these  objects  were  designated  as  planetary  candi¬ 


dates  KOI-701.01,  701.02,  and  701.03  in  the  Kepler 
2011  catalog  (7)  and  the  Kepler  2012  catalog  (8). 
KOI-701.04  and  701.05  were  subsequently  iden¬ 
tified  using  a  larger  data  sample  (9). 

Analysis  of  high-resolution  spectra  indicates 
that  Kepler-62  is  a  K2V  spectral  type  with  an 
estimated  mass  and  radius  (in  solar  units)  of 
0.69  ±  0.02  M0  and  0.63  ±  0.02  RQ  (9).  Exam¬ 
ination  of  the  sky  close  to  Kepler-62  showed  the 
presence  of  only  one  additional  star  that  con¬ 
tributed  as  much  as  1%  to  the  total  flux  (figs.  S3 
and  S4)  (9).  Warm-Spitzer  observations  (fig.  S9) 
and  the  analysis  of  centroid  motion  (table  SI) 
were  consistent  with  the  target  star  as  the  source 
of  the  transit  signals  (Fig.  1  and  fig.  SI).  We  com¬ 
puted  the  radius,  semimajor  axis,  and  radiative 
equilibrium  temperature  of  each  planet  (Table  1) 
on  the  basis  of  light  curve  modeling  given  the 
derived  stellar  parameters  (table  S3). 

The  masses  of  the  planets  could  not  be  di¬ 
rectly  determined  using  radial  velocity  (RV)  mea¬ 
surements  of  the  host  star  because  of  the  planets’ 
low  masses,  the  faintness  and  variability  of  the 
star,  and  the  level  of  instrument  noise.  In  the 
absence  of  a  detected  signal  in  the  RV  measure¬ 
ments  (9),  we  statistically  validated  the  plane¬ 
tary  nature  of  Kepler-62b  through  -62f  with  the 
BLENDER  procedure  (10-13)  by  comparing  the 
probability  of  eclipsing  binaries  and  other  false¬ 
positive  scenarios  to  bona  fide  transiting  planet 
signals  (14-18). 

To  systematically  explore  the  different  types 
of  false  positives  that  can  mimic  the  signals,  we 
generated  large  numbers  of  synthetic  light  curves 
that  blend  together  light  from  multiple  stars  and 
planets  over  a  wide  range  of  parameters,  and  then 
compared  each  blend  with  the  Kepler  photom¬ 
etry  (Fig.  2).  We  rejected  blends  that  resulted  in 
light  curves  inconsistent  with  the  observations. 
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